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ABSTRACT 

Context. Although there is currently a debate over the significance of the claimed large-scale anomalies in the Cosmic 
Microwave Background (CMB), their existence is not totally dismissed. In parallel to the debate over their statistical 
significance, recent work has also focussed on masks and secondary anisotropics as potential sources of these anomalies. 
Aims. In this work we investigate simultaneously the impact of the method used to account for masked regions as well 
as the impact of the integrated Sachs- Wolfe (ISW) effect, which is the large-scale secondary anisotropy most likely to 
affect the CMB anomalies. In this sense, our work is an update of both Francis & Peacock 2010 and Kim et al. 2012. 
Our aim is to identify trends in CMB data from different years and with different mask treatments. 
Methods. We reconstruct the ISW signal due to 2 Micron All-Sky Survey (2MASS) and NRAO VLA Sky Survey (NVSS) 
galaxies, effectively reconstructing the low-redshift ISW signal out to z ^ 1. We account for regions of missing data 
using the sparse inpainting technique of Abrial et al. 2008, Starck, Murtagh & Fadili 2010 and Starck, Fadili & Rassat 
2012. We test sparse inpainting of the CMB, Large Scale Structure and ISW and find that it constitutes a bias- free 
reconstruction method suitable to study large-scale statistical isotropy and the ISW effect. 

Results. We focus on three large-scale CMB anomalies: the low quadrupole, the quadrupole/octopole alignment, 
and the octopole planarity. After sparse inpainting, the low quadrupole becomes more anomalous, whilst the 
quadrupole/octopole alignment becomes less anomalous. However, after subtraction of the ISW signal, the trend 
amongst the CMB maps is that both the low quadrupole and the quadrupole/octopole alignment have reduced signif- 
icance, yet other hypotheses remain possible as well (e.g. exotic physics). Our results also suggest that both of these 
anomalies may be due to the quadrupole alone. While the octopole planarity significance is reduced after inpainting 
and after ISW subtraction, we do not find that it was very anomalous to start with. In the spirit of participating in 
reproducible research, we make all codes and resulting products which constitute main results of this paper public here: 
http : //www . cosmostat . org/anomaliesCMB . html. 



1. Introduction 

In recent years, cosmological observations (Larson et al. 
2011; Percival et al. 2007 a; Schrabback et al. 2010) 
have lead to the establishment of a standard cosmological 
model. This model assumes an inflationary scenario lead- 
ing to Gaussian features in the temperature anisotropics of 
the Cosmic Microwave Background (CMB). Since the ad- 
vent of the COsmic Background Explorer (COBE, Bennett 
et al. 1990) and Wilkinson Microwave Anisotropy Probe 
(WMAP, Spergel et al. 2003) several signatures of lack of 
statistical isotropy, or "anomalies", have been reported on 
large scales. 

On the largest scales, a low quadrupole was reported 
with COBE data (Hinshaw et al. 1996; Bond et al. 1998) 
and later confirmed with WMAP data (Spergel et al. 2003), 
suggesting that it was not due to Galactic emissions. The 
octopole presented an unusual planarity and a correlation 
with the quadrupole (Tegmark et al. 2003; de Oliveira- 
Costa et al. 2004; Slosar & Seljak 2004; Copi et al. 2010). 
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Other anomalies include a north/south power asymmetry 
(Eriksen et al. 2004; Bernui et al. 2006), an anomalous "cold 
spot" in the CMB (Vielva et al. 2004; Cruz et al. 2005; 
Cruz et al. 2006), alignments of other large-scale multi- 
poles (Schwarz et al. 2004; Copi et al. 2006), the so-called 
Axis of Evil (Land & Magueijo 2005a), and other violations 
of statistical isotropy (Hajian & Souradeep 2003; Land & 
Magueijo 2005b). 

If confirmed, these anomalies could provide a new win- 
dow into exotic early- Universe physics. However, there is 
much debate about the possible causes of these claimed 
anomalies. The statistics used to detect them are often 
subtle, and given the large cosmic variance on the scales 
considered, the anomalies could be due to a simple statisti- 
cal fluke or in fact not be anomalous (Bennett et al. 2011; 
Efstathiou et al. 2010; Gold et al. 2011; Bennett et al. 2012), 
depending on how the significance is measured. 

The anomalies could also have a low-redshift cosmologi- 
cal origin. Since statistical isotropy is predicted for the early 
Universe, analyses should focus on the primordial CMB, i.e. 
one from which secondary low-redshift cosmological signals 
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have been removed. This should be done whether or not 
one believes the significance of the reported anomalies in 
the CMB. 

Peiris & Smith (2010) investigated the kinetic Sunyaev- 
Zeldovich (kSZ) effect, a known secondary anisotropy, and 
found it was unlikely that the kSZ effect was the origin 
of these anomalies. On the largest scales, the only sec- 
ondary anisotropy is the integrated Sachs- Wolfe (ISW) ef- 
fect, which is correlated with foreground large-scale struc- 
ture. A first study by Rudnick et al. (2007) detected a cold 
spot in the NRAO VLA Sky Survey (NVSS), similar to that 
found in WMAP data, suggesting that the WMAP cold 
spot could be due to the late ISW effect. However, Smith 
& Huterer (2010) showed that the cold spot in the NVSS 
was no longer significant when systematics were taken into 
account, and Granett et al. (2009, 2010) rule out that the 
cold spot be due to a supervoid. Recently, Yershov et al. 
(2012) found a correlation between the spatial distribution 
of supernovae and the CMB, which they claim could be due 
to dust contribution, ionised gas, or the ISW effect. 

Francis & Peacock (2010) (hereafter FPIO) investigated 
the impact of the ISW field due to 2 Micron All-Sky 
Survey (2MASS) galaxies on several CMB anomalies: the 
low quadrupole, the quadrupole/octopole alignment, the 
planarity of the octopole, the north/south asymmetry and 
the cold spot. They found that the first two of these anoma- 
lies were reduced in significance after the removal of the 
foreground signal, but that the octopole planarity was actu- 
ally increased. The north/south asymmetry was somewhat 
reduced and the cold spot remained anomalous. 

The study of these large-scale anomalies is extremely 
sensitive to the treatment of the Galactic mask to ac- 
count for foreground removals (Bernui et al. 2006; Slosar 
& Seljak 2004; Bennett et al. 2012) in both the CMB and 
galaxy survey data used to reconstruct the ISW field, hence 
a statistically un-biased method for dealing with missing 
data should be used. Regarding the CMB, Bennett et al. 
(2012) found that the limiting factor in measuring the 
quadrupole/octopole alignment was foreground removal. 
Regarding the LSS side, FPIO reconstructed missing data 
by first filling the masked region with a random Poisson 
sampling of galaxies with the same average number den- 
sity as outside the mask and then applying a Wiener filter, 
which is optimal for Gaussian data. Kim et al. (2012) in- 
vestigated methods to reconstruct CMB data in masked 
regions using a Gaussian-constrained harmonic inpainting 
method. They found that the quadrupole/octopole align- 
ment in CMB maps was increased after inpainting treat- 
ment. Both methods assume Gaussianity of the underlying 
maps, which is a limitation in the context of search for sta- 
tistical anisotropy. The filling of the masked region with 
random galaxies may also cause an artificial quadrupole 
since the Galactic mask has the shape of a quadrupole. 

In this paper, we investigate the use of sparse inpaint- 
ing (see Abrial et al. (2008) and Starck, Murtagh, & Fadih 
(2010), hereafter A08-SMF10) on CMB (WMAP) and LSS 
maps (2MASS and NVSS) in the context of probing large- 
scale anomalies in the CMB. The sparse inpainting does not 
assume that the true underlying map is Gaussian, or sta- 
tistically isotropic. We reconstruct a full-sky tomographic 
ISW map due to 2MASS and NVSS galaxies. We test 
whether removing the ISW signal due to these LSS maps 
affects the observed large-scale anomalies in the CMB by 
focussing on the same three large-scale anomalies as FPIO. 



This paper updates the work presented in Francis & 
Peacock (2010) and Kim et al. (2012), with the following 
differences: 

1. Tomography: We reconstruct the ISW signal due to both 
2MASS (Jarrett et al. 2000) and NVSS (Condon et al. 
1998) data. The data sets are described in Section 3 

2. Trends: We look for trends by considering 11 different 
WMAP data sets, described in Section 3 

3. Sparse inpainting: Missing data in CMB and LSS maps 
are reconstructed using the sparse inpainting technique 
of A08-SMF10 and Starck, Fadih, & Rassat (2012), de- 
scribed in Section 4, which does not assume that the 
data is Gaussian or statistically isotropic 

4. Tests for biases: We test that the sparse inpainting 
method does not introduce biases in statistical isotropy 
tests nor creates a spurious ISW signal in Section 4.1 

In Section 2 we present how it is possible to estimate 
the primordial CMB on large scales with knowledge of 
foreground LSS maps by reconstructing the ISW field. In 
Section 5, we investigate the effect of subtracting the ISW 
field due to 2MASS and NVSS galaxies on three previ- 
ously reported large-scale anomalies: the low quadrupole, 
the quadrupole/octopole alignment and the octopole pla- 
narity. In Section 6, we present a discussion of our results. 

2. Estimating the Large-Scale Primordial CMB 

2.1. The large-scale primordial CMB after subtraction of the 
integrated Sachs-Wolfe effect 

The observed temperature anisotropics in the CMB, Sqbs^ 
can be described as the sum of several components: 

^^OBS = ^^prim + ^ISW + ^other + A/" , (1) 

where ^prim are the primordial temperature anisotropics, 
^isw total secondary temperature anisotropics due 

to the ISW effect, (Mother are other secondary anisotropics 
(e.g., Sunyaev-Zeldovich effect), and Af is noise. Here, we 
assume that the observed data are clean of any foregrounds 
and that any regions requiring masking due to contami- 
nated data will be corrected for during the sparse inpainting 
phase. On large scales, the ISW signal is the only secondary 
anisotropy and the anisotropics are cosmic variance limited, 
so that the last two terms can be ignored: 

Sobs ^ ^^pdm + for large scales. (2) 

The ISW effect arises in universes where the cosmic poten- 
tial decays at late times, as is the case with dark energy, 
open curvature, or possibly some modified gravities. The 
temperature anisotropics due to the ISW effect are given 
by: 

<5i'sw=-2/ $'((%-7?)n,r,)d7?, (3) 

where 77 is the conformal time, defined by dr] = and 
7^0 and r]L represent the conformal times today and at the 
surface of last scattering respectively. The unit vector n 
is along the line of sight and the gravitational potential 
<l>(x, T]) depends on position and time. The integral depends 
on the rate of change of the potential = d^/dr]. The 
potential field can be related to the matter field, of which 
a galaxy map is assumed to be a tracer. 
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Table 1. : List of the 11 teraperature anisotropy maps we use to probe the anomalies in the primoridal CMB, including 
year of WMAP data they correspond to and a description of the mask treatment. 'TOH' corresponds to the treatment 
in Tegmark et al. (2003), TLC to the Internal Linear Combination method (see Hinshaw et al. 2007; Bennett et al. 2012; 
Gold et al. 2009, 2011), TLC- Wavelets' corresponds to the method in Delabrouille et al. (2009), and the sparse inpainting 
method is described in A08-SMF10 and uses the sparsity prior described in Starck, Fadili, & Rassat (2012). 



By subtracting the reconstructed ISW signal, (5isw, due 
to the matter field traced by a foreground galaxy survey, we 
can estimate the large-scale primordial CMB temperature 
anisotropy field 6prim by: 



-'prim 



^^OBS — (^isw, for large scales. 



(4) 



where this relation tends to equality (on large scales) when 
the entire mass distribution of the Universe is used to re- 
construct the ISW field, i.e. when ^isw = ^isw • 

In this paper, we estimate (5prim by subtracting the ISW 
temperature anisotropy field as reconstructed in Section 4 
using 2MASS and NVSS data (see Section 3): 



^OBS 



r2MASS 
^ISW 



rNVSS 
^ISW 



(5) 



where the extra last term corresponds to the removal of 
the kinetic Doppler quadrupole (ignoring the monopole and 
dipole). The kinetic Doppler quadrupole contribution is 
given by: 

cos^O-^ , (6) 



'JkD,£ 



where is the angle between the position on the sky 
and the direction of motion creating the kinetic Doppler 
quadrupole (Copi et al. 2006). For our calculations, we 
take V = 370km5-^ towards (/,6) = (263.85°, 48.25°), 
and c is the speed of light, similarly to what is done in 
Francis & Peacock (2010). We make the kinetic Doppler 
quadrupole map in Healpix format as well as code to gener- 
ate it available for download here http : //www . cosmostat . 
org/ anomaliesCMB . html. 

2.2. The ISW temperature field from LSS maps 

The temperature ISW field can be reconstructed in spheri- 
cal harmonics, ^ from the LSS field g^rn (Boughn et al. 
1998; Cabre et al. 2007; Giannantonio et al. 2008): 



cISW _ ^ - - ^ 



(7) 



where gim represent the spherical harmonic coefficients of 
a galaxy overdensity field g{0^ </>), given by 



^(6>,0) = ^^£m>£m(6>,^), 



(8) 



where Yim{0^(j)) are the spherical harmonics. The spectra 
Cgg and CgT are the galaxy (g) and CMB (T) auto- and 
cross-correlations measured from the data or their theoret- 
ical values given by: 

C,T W = 47r6, j dk^Wg(k)WT(k), (9) 



where 



Wrik) 



Wg{k) = / dre(r)j,(kr)D(z), 



(10) 
(11) 



drj,(kr)H(z)D(z)(f-l), (12) 



e(r) 



(27r)3 



(13) 



(14) 



J drr^n(r) ' 

where we use the same notation as in Rassat et al. (2007) 
and have assumed a linear bias b{k, z) = hg and D{z) is the 
linear growth function. 

From Equation 7, we can check that: 

(C'^sL) = c,T{e). (15) 

and 



/ rISW r*ISW\ _ 
\^lm ^im / — 



C'gg(^) 



= Cisw(^). 



(16) 



Equation 7 and 16 show that the ISW temperature field is 
independent of the galaxy bias 6^, which relates the galaxy 
and matter (m) fluctuations by g^rn = ^g^m,imi and that it 
is not necessary to estimate the value of the galaxy bias in 
order to estimate the ISW temperature field. 

We note that the above method is suitable for photo- 
metric surveys where a projected spherical harmonic de- 
composition is satisfactory. For spectroscopic surveys, a full 
three-dimensional reconstruction can be done in spherical 
Fourier-Bessel decomposition, as in Shapiro et al. (2012) 
and Rassat & Refregier (2012) for the theoretical calcula- 
tions, and using a code such as 3DEX (Leistedt et al. 2012) 
for the measurements. 
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WMAP9 Masked 



WMAP9 Inpainted 





3MASS Masked 



3MASS Inpainted 




Fig. 1. WMAP ILC 9th year CMB data (top in mK), 2MASS data (middle), and NVSS data (bottom) with mask (left) 
and after sparse inpainting (right). Data in the right column is presented up to ^ = 10 for all maps. The original 2MASS 
data (middle left) is plotted with a maximum overdensity value of 3 to increase contrast in the map. 



3. Data 

We are interested in reconstructing the temperature ISW 
fields due to 2MASS and NVSS data. We describe the CMB 
and LSS data we use in the following section, while the 
reconstruction is described in Section 4. 



3.1. Cosmic Microwave Background Data 

For the CMB, we investigate data from several years of 
WMAP, including independent treatments of masked data. 
We do this to identify trends in the observed large-scale 
anomalies. We use 11 maps in total: : the Tegmark et al. 
(2003) reduced-foreground CMB map (TOH); the Internal 
Linear Combination Map (ILC) WMAP data from the 3rd 
year (Hinshaw et al. 2007) \ 5th year (Gold et al. 2009) 2, 



^ http: //lambda. gsfc .nasa.gov/product/map/dr2/ilc_ 
map_get . cf m 

^ http : //lambda . gsf c . nasa . gov/ product /map/ dr3/ ilc_ 
map_get . cf m 



7th year (Gold et al. 2011)^, and 9th year (Bennett et al. 
2012) ^ as well as sparsely inpainted versions of the ILC and 
TOH maps (see Section 4). We also include the sparsely 
inpainted WMAP 5th year ILC map by Delabrouille et al. 
(2009), which was reconstructed using a wavelet technique. 
We summarise the CMB maps used in Table 1. The ILC 
9th year temperature overdensity map is shown with its 
corresponding mask ^ in the left-hand side of Figure 1 (top). 

3.2. The 2 Micron All Sky Survey (2MASS) Data 

As a tracer of the low-redshift matter distribution, we use 
the publicly available 2MASS full-sky extended source cat- 



http : //lambda . gsf c . nasa . gov/ product /map/ dr4/ ilc_ 
map_get . cf m 

^ http: //lambda. gsfc .nasa.gov/product/map/dr5/ilc_ 
map_get . cf m 

^ We use the 7th year temperature analysis mask for 
nside — 512 (http: //0-lambda. gsf c .nasa.gov. iii-server . 
ualr . edu/product/map/dr4/masks_get . cf m) 
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alogue (XSC) selected in the near IR (Jarrett et al. 2000), 
which has median redshift z ~ 0.07. The near-IR selection 
means galaxies are relatively well observed, even in the re- 
gion of the Galactic plane. As a result, the 2MASS survey 
has a large sky coverage, which is ideal for studying the 
ISW effect, with j^^^^^ = 0.69. For exact details of our 
selection criteria and the mask, see Rassat et al. (2007) 
and Section 6.2 of Dupe et al. (2010). The 2MASS over- 
density map is shown with its corresponding mask on the 
left-hand side of Figure 1 (middle). 

3.3. The NRAO VLA Sky Survey (NVSS) Data 

As a tracer of the intermediate redshift matter distribu- 
tion, we use the publicly available NRAO VLA Sky Survey 
(NVSS), which has a median redshift z - 1.4. The NVSS 
survey is a 1.4GHz (radio) continuum survey covering the 
sky north of S = —37°, in which most sources away from 
the Galactic plane are of extragalactic origin (Smith, Zahn, 
& Dore 2007). We select only sources with flux larger than 
2.5 mJy. For the mask we exclude the Galactic region with 
\b\ < 10° and also mask out a 0.6° radius around bright 
sources with flux larger than 2.5 Jy, as in Ho et al. (2008). 
The resulting mask gives us an effective sky fraction of 
/NVSS ^0.66. 

In order to account for the problem of declination- 
dependent density, which NVSS data suffers from, we apply 
the same method as in Schiavon et al. (2012) and Vielva 
et al. (2006): the NVSS overdensity is calculated separately 
for nine declination bins with Asin((5) = 0.1. 

Alternatives to the declination problem have been 
proposed, e.g. by considering only objects above 5mJy 
(Barreiro et al. 2012) or including templates to project out 
the declination-striping modes (Ho et al. 2008). However, 
none of these methods satisfactorily solves a problem 
of missing power on large scales for NVSS (Hernandez- 
Monteagudo 2010). The NVSS overdensity map is shown 
with its corresponding mask on the left-hand side of Figure 
1 (bottom). 

4. Reconstruction of Full-Sky Maps 

Large-scale modes in the CMB and LSS fields are very sen- 
sitive to large coherent regions of missing data, such as 
those which are due to a Galactic mask. Therefore it is 
crucial to use a reconstruction method which does not in- 
troduce biases in the reconstructed field. To account for re- 
gions of missing data in CMB and LSS maps, we use sparse 
inpainting (A08-SMF10) to reconstruct regions of missing 
data (see also Appendix A in Dupe et al. 2010). The sparse 
inpainting approach is powerful as the only assumption it 
makes on the underlying field is sparse representation. In 
this particular case, the assumption is therefore that the 
CMB and LSS signals are sparsely represented in spherical 
harmonic space, i.e. only a few a^rnS (not C{i)s) are re- 
quired to describe the data. This is easily verified directly 
from the data and is described in more detail in Starck et al. 
(2012) and Dupe et al. (2010). 

Sparse inpainting does not assume statistical isotropy or 
Gaussianity of the underlying field, unlike other methods 
such as Wiener filtering (used by Francis & Peacock 2010) 
or constrained Gaussian realisations (as used by Kim et al. 
2012). 



In section 4.1, we test whether the sparse inpainting re- 
construction might produce biases in measurements of sta- 
tistical isotropy or a spurious ISW signal. In section 4.2, we 
describe the details of the sparse inpainting reconstruction 
for the LSS and CMB maps and in section 4.3, we describe 
how we use these to reconstruct full-sky ISW maps. 

4.1. Sparse Inpainting as a Bias- Free Reconstruction Method 

In order to use the reconstructed maps to test large-scale 
anomalies in the CMB, we must first be confident that there 
are no spurious large-scale anomalies due to our reconstruc- 
tion technique. The full details of the tests are given in 
Appendix A, and we highlight the main conclusions here. 

We first test whether sparse inpainting affects the three 
anomalies we are testing in the CMB data (see Section 5 for 
the anomaly descriptions). We test this by considering two 
sets of Gaussian random field simulations of CMB data. 
The first set of simulations has a low quadrupole as the- 
oretical input, and the second set has a WMAP7 best-fit 
theoretical quadrupole (see Equation 17). We test for the 
three anomalies for both sets of simulations before and af- 
ter inpainting (Tables A. 2, A.l). We also apply the same 
tests to Gaussian random field simulations of 2MASS and 
NVSS galaxy fields, using theoretical values of their spectra 
(including the galaxy bias. Tables A. 3 and A. 4). 

We find that the sparse inpainting method applied to 
statistically isotropic CMB and LSS simulated data does 
not introduce any significant biases for any of the three 
tests of statistical isotropy. After inpainting, the quadrupole 
is slightly lower but only by 3%, whether considering the 
simulations with low quadrupole or theoretical quadrupole 
on input. The test of quadrupole/octopole alignment can be 
somewhat altered due to sparse inpainting (see the standard 
deviation on the bias). We show here that if the CMB and 
LSS fields are statistically isotropic, then sparse inpainting 
does not alter this. 

Finally, we want to test whether sparse inpainting will 
alter in any way the reconstructed ISW temperature sig- 
nal, either by creating a spurious ISW signal or biasing the 
tests of statistical isotropy. To do this, we first consider 
simulated NVSS maps with uncorrelated CMB maps (1st 
line of Table A. 5). We test for cross-correlations before and 
after inpainting of the individual NVSS and CMB maps 
(which have different masks). Since the maps are statis- 
tically uncorrelated, we expect a null cross-correlation on 
average, which we find is the case before and after applying 
sparse inpainting. This shows that sparse inpainting does 
not introduce a spurious ISW or cross-correlation signal, 
confirming what we had found in Dupe et al. (2010). 

We also test the quadrupole/octopole alignment and 
octopole planarity of the reconstructed ISW signal from 
NVSS galaxies before and after inpainting (again where the 
NVSS and CMB maps are inpainting separately with dif- 
ferent masks, see Table A. 5) and find that the tests are 
unbiased after sparse inpainting. 

We note that it would, of course, be interesting to know 
whether sparse inpainting is also unbiased in the case of 
statistically anisotropic fields, for example, the case where 
the CMB really does have anomalies before sparse inpaint- 
ing is applied. However to do this, one would first have to 
decide on the models causing the anomalies (exotic physics, 
foregrounds, etc. . . ) in order to test realistic statistical 
anisotropy. In this paper, we focus only on finding out if 
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sparse inpainting creates spurious statistically anisotropic 
signatures in data where the true distribution is statistically 
isotropic. 

4.2. Full-sky CMB and LSS maps 

We describe here the details for the sparse inpainting recon- 
struction for the CMB, 2MASS and NVSS data. Our goal is 
to reconstruct the quadrupole and octopole only {£ = 2, 3), 
since we are only interested in anomalies on the largest 
scales. However, due to mode correlations induced by the 
masked data, we must consider multipoles larger than £ > 3 
for the reconstruction. 

For the 11 CMB maps we reconstruct the harmonic 
coefficients up to = 64 using nside = 64 using the spar- 
sity prior as described in Starck, Fadili, & Rassat (2012), 
i.e. using the following command line in the open-source 
sparse inpainting package I SAP software^: 

> aim = cinb_lowl_alin_inpainting(inap. Mask, 
linax=linax, niter=100, InpMap=result) . 

For 2MASS data, we reconstruct the harmonic coeffi- 
cients up to ^ = 64 using nside = 64 with the following 
command line in I SAP software^: 

> aim = cmb_lowl_alm_inpainting(t , 
Mask , lmax=lmax , niter=500, InpMap=result , 
/galaxy) , 

where the keyword /galaxy was optimised for point 
source catalogues like galaxy surveys to take into account 
that the field does not have a zero mean. 

For NVSS data, we reconstruct the harmonics up to 
^ = 64 using nside = 128 and the same command line and 
options as for 2MASS. The WMAP 9th year ILC, 2MASS 
and NVSS maps are shown before (left) and after inpainting 
(right) in Figure 1. The ten other CMB maps are also re- 
constructed in the same manner, but not plotted in Figure 
1. 

In Figure 2, we plot the measured quadrupole and oc- 
topole from all CMB maps listed in Table 1 before (black) 
and after (red) sparse inpainting. Apart from the TOH Wl 
map, the values of the quadrupole and octopole are similar 
across different WMAP years. After inpainting, both the 
quadrupole and octopole are lowered for all maps, with the 
TOH Wl map again having slightly different values than 
the other maps for the quadrupole. 

The simulations in Appendix A show that inpainting 
on W7-like maps (i.e. with low quadrupole, Ci=2 = 250.6 ± 
161.0 fiK'^) introduces only a slight bias (7.2, i.e. about 
a 3% drop) on the quadrupole, so that we do not expect 
the drop in quadrupole value shown in Figure 1 to be an 
artifact of inpainting. 

Regarding the LSS data, we find that after reconstruc- 
tion of the NVSS map the theoretical power spectrum as 
predicted from the redshift distribution of radio sources 
given by Ho et al. (2008) does not correspond well to the 
measured power spectrum, especially at higher multipoles 
(using their value of the galaxy bias 6nvss = 1.98), as al- 
ready noted by Hernandez-Monteagudo (2010). However, 
since we are focussing on very large angles only {i = 2, 3), 
it is difficult to assess whether the theoretical spectrum 
needs to be revisited. For this work, we use the N{z) pro- 
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Fig. 2. Quadrupole and octopole values {fiK'^) for the dif- 
ferent CMB maps listed in Table 1. Points in black are with 
no extra treatment than that given in the literature, and 
points in red are after application of sparse inpainting. 




Fig. 3. Redshift distributions for 2MASS (see Afshordi 
et al. 2004; Rassat et al. 2007; Dupe et al. 2010) and NVSS 
(from Ho et al. 2008). 



vided from Ho et al. (2008), which is only used in the NVSS 
simulations for testing sparse inpainting (Appendix A) and 
for the prediction of the value of quadrupole and octopole 
of the ISW due to NVSS data (Table 2 in Section 4.3). The 
actual reconstructed ISW map uses the measured auto- and 
cross- spectra (see the following Section 4.3) and is therefore 
independent of theory. 

4.3. Full-sky ISW temperature maps 

We reconstruct the full-sky ISW maps due to 2MASS and 
NVSS data from their inpainted maps (Figure 1), using 
Equation 7. Since the redshift distribution of 2MASS and 
NVSS has little overlap (see Figure 3), we calculate the 
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total ISW signal simply by summing the ISW maps recon- 
structed from each survey individually. 
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Table 2. Amplitude of the ISW temperature quadrupole 
and octopole due to the 2MASS galaxy survey {left: as 
measured from data using inpainted WMAP 9 data for 
the CMB, and right: from theoretical calculations) and a 
Euclid- type survey (using theoretical spectra). The stan- 
dard deviations are calculated assuming f^^^^^ = 0.69, 



/,^jfss ^ 0.66 and /™ = 0.48. 
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Fig. 4. Quadrupole and octopole values (/ii^^) for the dif- 
ferent reconstructed ISW maps. The ISW maps are recon- 
structed from sparsely inpainted CMB and LSS maps using 
Equation 7, where the auto- and cross-correlations are mea- 
sured directly from the data. Small variations in the phases 
of different renditions of CMB maps can lead to negligi- 
ble differences in the auto-correlation (see Figure 2) and to 
larger differences in a cross-correlation with a second map. 
This in turn can result in larger variations in the ampli- 
tude of the ISW signal, explaining why the reconstructed 
ISW maps vary slightly more between CMB maps than the 
observed temperature power does. 

We estimate the cross- and auto-correlations in 
Equation 7 using data only with inpainted maps for both 



the LSS and CMB data. This has the advantage of assum- 
ing only the form of the ISW signal, but not its presence, 
i.e. in the absence of a correlation between the CMB and 
the LSS maps, the reconstructed ISW map would be zero. 
As mentioned in Section 2.2, no numerical knowledge of 
the linear bias is necessary in this case. This method also 
means there is an ISW map for each CMB map considered. 

Quadrupole and octopole values for the reconstructed 
maps using WMAP 9 ILC inpainted map are reported in 
Table 2 and plotted for each CMB map in Figure 4. 

The theoretical values in Table 2 are calculated employ- 
ing the N{z) used in Afshordi et al. (2004) and Rassat et al. 
(2007) for 2MASS data and Ho et al. (2008) for NVSS data. 
They also require the assumption of a fiducial cosmologi- 
cal model, for which we assume a "vanilla" model, not a 
specific best-fit value, given that we are considering several 
CMB maps. The vanilla model we chose has the following 
cosmological parameters: = 0.25, ft^ = 0.045, I^de = 
0.75, k;o = — = 0,ns = l,/i = 0.7,(78 = 0.80, i.e. a 
universe which can have curvature and does not have mas- 
sive neutrinos. 

The theoretical values in Table 2 are calculated us- 
ing the exact formulae (not the Limber approximation). 
The error bars due to cosmic variance are calculated using 
_^2MASS^0.69, and /NVSS ^0.66. 

For comparison, we also indicate the expected value 
of the ISW map from Euclid data (Laureijs et al. 2011; 
Refregier et al. 2010). We predict this by considering a 
galaxy survey with median redshift z = 0.80 and a Smail- 
type redshift distribution with a = 2,/3 = 1.5 (see for e.g. 
Equation 19 in Kirk et al. 2011) and /sky = 0.48. We calcu- 
late the expected ISW signal assuming a single large red- 
shift bin. 

Figure 4 shows the quadrupole and octopole power for 
the six reconstructed ISW maps (we only consider inpainted 
maps for the ISW reconstruction). The ISW quadrupole 
varies from 96.08 - 117.13/17^^ (except for TOH, which 
returns a significantly larger quadrupole of 149.07/ii^^, 
as shown in Figure 4), while the octopole varies from 
37.70 - 47.13/ii^^, depending on the WMAP inpainted 
data used. Since each map is estimated from the cross- 
correlation of the reconstructed LSS maps with the cor- 
responding reconstructed CMB map, this may explain the 
larger variation in Figure 4 than in Figure 2. This is because 
a small phase change may not change the auto-correlation 
(as in Figure 2), but can lead to larger variations in cross- 
correlations with another map. This in turn can lead to 
larger variations in the amplitude of the ISW signal. 

In our analysis, we chose to use all six reconstructed 
ISW maps, because each reconstructed ISW map is self- 
consistently produced, i.e. without any assumption about 
the theoretical amplitude of the ISW signal. In addition, 
we are interested in identifying possible trends in the way 
subtraction of the ISW signal affects the CMB anomalies, 
which will provide a stronger case for any conclusions. 

Finally, Figure 5 shows the reconstructed ISW tempera- 
ture maps due to 2MASS galaxies (top) and NVSS galaxies 
(middle) where the CMB map is the W7 inpainted data. 
Only the quadrupole (left) and octopole (right) are shown, 
since we are interested in probing large-scale anomalies at 
^ = 2 and ^ = 3. 

The quadrupole and octopole of the ISW map due to 
2MASS galaxies can be compared with those obtained by 
Francis & Peacock (2010) (top of their Figure 3). The main 
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difference is that our reconstructed quadrupole is much less 
planar. This may be due to the difference in reconstruction 
of missing data. Francis & Peacock (2010) reconstructed 
missing data by first filling the masked region with a ran- 
dom Poisson sampling of galaxies with the same average 
number density as outside the mask. However, this can 
cause an artificial quadrupole, since the Galactic mask has 
the shape of a quadrupole. 

5. Probing Anomalies in the Observed and 
Primordial CMB 

The main goal of this paper is to test whether reported 
large-scale anomalies in the observed CMB are still present 
in the primordial CMB, which we estimate by subtract- 
ing the ISW signal. Since the estimation of the ISW signal 
requires sparse inpainting of all CMB and LSS maps to ac- 
count for regions of missing data, we can also test whether 
reported large-scale anomalies persist in the observed CMB 
after inpainting. 

We investigate three large-scale anomalies in the CMB, 
which are the same as those investigated in Francis & 
Peacock (2010). In Section 5.1, we consider the impact on 
the low quadupole power, in Section 5.2 the alignment of 
the quadrupole and octopole, and in Section 5.3 the pla- 
narity of the octopole. 

5.1. Low Quadrupole Power 

The low power of the quadrupole was first reported for 
COBE data (Bennett et al. 1992; Bond et al. 1998; Hinshaw 
et al. 1996) and subsequently observed in WMAP data 
(Spergel et al. 2003). Depending on how the significance 
of the low quadrupole is measured, recent WMAP papers 
now argue that there is no anomaly in the amplitude of 
the quadrupole (Bennett et al. 2011, 2012). However, we 
are still interested in studying how the amplitude of the 
quadrupole behaves after sparse inpainting and after sub- 
traction of the reconstructed ISW field. 

The values of the quadrupole before subtraction of the 
ISW field are reported in the top part of Table 3 [labelled 
'1)' and '2)'] before and after sparse inpainting. The third 
column of Table 3 gives the probability ^ for a random 
variable with 5 degrees of freedom to take a value less than 
or equal to the WMAP9 (Bennett et al. 2012) expected 
theoretical (TH) value: 

Cw^/=2 = 1161.33 fiK^. (17) 

The inpainted maps, which were shown not to introduce 
a bias in the quadrupole reconstruction (see Section 4 and 
Appendix A), have lower quadrupole values than the input 
maps. This means that the application of inpainting for the 
CMB reconstruction actually increases the significance of 
the low quadrupole anomaly for all considered CMB maps. 

The third part of Table 3 [labelled '3)'] shows the mea- 
sured quadrupole power and probability after subtraction 
of the ISW reconstructed map due to 2MASS and NVSS 

^ This is done using the idl routine 
chisqr_pdf (quad*df /theory , df), where quad is 
the observed value of the quadrupole, theory the 
theoretical value given by Equation 17, and df the 
number of degrees of freedom for a quadrupole, i.e. 5. 
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Table 3. Quadrupole {£ = 2) power and corresponding 
probability of the quadrupole power being so low. 1): For 
the best fit of WMAP 7th year and 9th year data. 2): For 
11 different CMB maps. 3): After subtraction of the re- 
constructed ISW signal due to 2MASS and NVSS galaxies. 
4): Results from FPIO (Francis & Peacock 2010, bottom 
row). Probabilities in this work are calculated using the ex- 
pected theoretical value given in the last column taken from 
WMAP 9 best-fit results, whereas FPIO use best-fit results 
from WMAP 3. 



galaxies. The last column shows the new expected theoret- 
ical value of the CMB quadrupole, given by: 

/^TH _ r^TYL ^TH ^TH /-i q\ 

^£=2 — ^W9,£^2 ~ '-^ISW,2MASS,£=2 ~ ^ISW ,NVSS,£=2^ \^^) 

The values for C'jg^ 2mass/=2 ^isw,nvss,-^=2 
given by the values in Table 2 for the combined NVSS + 
2MASS case. These values are compared with those calcu- 
lated in FPIO [labelled '4)'] of Table 3. 

After ISW subtraction, the quadrupole power in- 
creases for ah CMB maps (from 102.1 - 124.7/iK^ to 
315.5 — 371.8/ii^^, omitting TOH, which has a much larger 
quadrupole value than all the other maps after inpainting), 
but not as significantly as it did in FPIO (from 250.6 to 

Subtracting the ISW from the observed CMB has two 
effects: the first is to reduce the theoretical quadrupole (this 
is always the case, for any non-zero ISW field). The second 
effect is that the measured quadrupole is larger (by chance). 
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Fig. 5. Quadrupole (left) and octopole (right) of the reconstructed ISW map for 2MASS (top) and NVSS (middle) and 
both NVSS and 2MASS together (bottom). The quadrupole and octopole of the ISW map due to 2MASS galaxies can 
be compared with that obtained by Francis & Peacock (2010) (top of their Figure 3). Maps are shown in units of jiiK. 



Both these effects reduce the anomaly independently. 
They explain why, even though the quadrupole after ISW 
subtraction is lower in our case than in FPIO, the prob- 
ability of having such a low quadrupole is still larger for 
some CMB maps. After subtraction the significance of the 
quadrupole anomaly is raised in the range 14.7 — 33.2%, 
instead of 0.58 - 4.4% before ISW subtraction. 



We note that the WMAP9 team uses a different ap- 
proach to measuring the significance of the low quadrupole 
(Bennett et al. 2012). In any case, our study shows that af- 
ter ISW subtraction the quadrupole cannot be considered 
low, whichever statistical approach is used. 



5.2. Quadrupole/Octopole Alignment 

It was first noted by de Oliveira-Costa et al. (2004) that 
not only did the CMB quadrupole and octopole both ap- 
pear planar (i.e. the dominant mode was m = ^), they also 
seemed aligned along a similar axis. For a Gaussian random 
field there is no reason for the preferred axes of two differ- 
ent multipoles to be correlated. For each multipole, this 
preferred axis can be quantified by maximising the quan- 
tity 



Kn) = V'm^|a£^(n)p, 



(19) 



where a£rn(n) corresponds to the a£m coefficients of the ro- 
tated CMB temperature anisotropy map, where n corre- 
sponds to the new z-axis, and is best calculated with at 
least nside = 512. The axes n2 and 113 correspond to the 
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Fig. 6. Quadrupole (left) and octopole (right) of WMAP 9 data before inpainting (top), after sparse inpainting (middle), 
and after subtraction of the reconstructed ISW signal due to 2MASS and NVSS galaxies (bottom). 



preferred axes (i.e. where g^(n) is largest) for i = 2 and 
£ = 3 respectively. We note that this way of measuring the 
quadrupole/octopole alignment differs from that used by 
Bennett et al. (2011, 2012), who maximize the quantity: 



(20) 



Bennett et al. (2012) found that the limiting fac- 
tor in measuring the quadrupole/octopole alignment was 
foreground removal and treatment of the Galactic mask. 
This motivated us to investigate how the measured 
quadrupole/octopole behaves after sparse inpainting. 

In the top part of Table 4, we present the dot prod- 
uct of the preferred axes n2, ns for ^ = 2, 3 for each CMB 
map before and after sparse inpainting and before and after 
ISW subtraction, as well as their corresponding separation 
(note the maximum separation is 90° since the axes are not 



vectors). We also note the probability of having such a sepa- 
ration in a Gaussian random field, which we estimate using 
the fact that the dot product of the two vectors has a uni- 
form distribution (de Oliveira-Costa et al. 2004). Finally, 
we compare this with the effect found by FPIO using their 
2MASS-ISW reconstruction subtracted from the ILC W3 
map (last line of Table 4). 

For non-inpainted maps, the separation varies from 2.3- 
9.7°, i.e. probabilities ranging from 0.08% to 1.4%. The 
same CMB maps, when inpainted, are less anomalous with 
separations ranging from 10.6-17.9°, corresponding to prob- 
abilities of 1.7-4.9% 

When the ISW field due to 2MASS and NVSS galax- 
ies is removed, FPIO noted that the alignment was signif- 
icantly reduced (see the last line in Table 4 or their Table 
2 for more details). Using our reconstruction, we also find 
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Table 4. Scalar product of the preferred axes of the 
quadrupole and octopole (n2 • ^3), its corresponding sep- 
aration (°), and the probability (%) of having such a low 
separation. Note: the theoretically allowed range is [0°— 90°] 
since the axes are not vectors. The results are compared 
with those of FPIO (Francis & Peacock 2010, bottom row). 



5.3. Planarity of the Octopole 

The third anomaly we investigate is the reported planarity 
of the CMB octopole, i.e. the fact that the phase m = ^ is 
preferred, which was first noted by de Oliveira-Costa et al. 
(2004). To quantify if the octopole is planar, de Oliveira- 
Costa et al. (2004) suggested measuring the quantity: 



|a3-3(n)P + |a33(n)|^ 



|tt3m(n) 



(21) 



This quantity represents the ratio of the octopole power 
which is contained in the mode m = ^ = 3, i.e. it is a test 
of planarity. For a Gaussian random field, the distribution 
of power amongst modes should be random, and so there 
is no reason that the value of t should be close to 1. 

In Table 5, we report values of the 't' statistic for CMB 
maps, before and after inpainting, and before and after ISW 
subtraction. The 't' value is calculated using nside = 128, 
and its corresponding probability is estimated using 1000 
simulated Gaussian random fields with the same power as 
the map considered. 

Before inpainting, the probability of having such planar 
octopoles ranges from 9.60-14.3%. de Oliveira-Costa et al. 
(2004) reported t = 0.94 and a probability around 7% us- 
ing the TOH map, so we first note that later WMAP data 
are less anomalous than the first year. After inpainting, oc- 
topole planarity is reduced even further, with probabilities 
ranging from 13.6 — 26.1%. After ISW subtraction, FPIO 
had found that the octopole became even more planar. We 
find the opposite, i.e. that after ISW subtraction, there is 
no evidence for any octopole planarity. However, we keep 
in mind that even before ISW subtraction the significance 
of the octopole planarity was not particularly significant. 



the alignment is reduced, though the significance is not as 
large as in FPIO. We note that we are removing the ISW 
from both 2MASS and NVSS galaxies, whereas FPIO only 
removed the ISW field due to 2MASS galaxies. Therefore, 
the general trend is that removing the ISW field reduces 
the anomalous alignment of the quadrupole and octopole. 

Figure 6 shows the quadrupole (left) and octopole 
(right) for the specific case of WMAP9 data, before in- 
painting (top), after inpainting (middle) and after ISW sub- 
traction (bottom). Inpainting reduces the power in both 
the WMAP9 quadrupole and octopole, as we had seen 
in Figure 2. Before inpainting, the preferred axes are in 
the direction of {l,b) = (-124.1,58.1) and (-122.3,62.9) 
for £ = 2,3 respectively. After inpainting, these become 
{l,b) = (-134.5,54.8) and (-110.8,57.7). After inpainting 
and subtraction of the ISW signal due to 2MASS and NVSS 
galaxies, the preferred axes become (/,6) = (—16.3,75.1) 
and (—86.5, 70.0) for ^ = 2,3 respectively. In other words, 
the quadrupole axis is the one which changes the most and 
which causes the anomaly to decrease since this anomaly 
measures the correlation between the preferred direction of 
both the quadrupole and the octopole. We find this trend is 
true for all maps considered, i.e. that it is the change to the 
quadrupole shape which is the main cause of the anomaly 
decrease, whether after sparse inpainting only or after ISW 
subtraction. 



6. Discussion 

Although there is currently a debate over the significance 
of the claimed large-scale anomalies in the CMB (see e.g., 
Bennett et al. 2011), their existence is not totally dis- 
missed. In parallel to the debate on the statistical signif- 
icance of the anomalies, recent work has also focussed on 
the impact of the reconstruction method to account for 
masked regions of the sky (e.g., Bernui et al. 2006; Slosar & 
Seljak 2004; Kim et al. 2012). Moreover, in some cases the 
Galactic mask is a limiting factor in the study of anomalies 
(Bennett et al. 2012). Studies have also focussed on low- 
redshift cosmology and astrophysics as potential sources of 
contamination (Rudnick et al. 2007; Peiris & Smith 2010; 
Francis & Peacock 2010; Yershov et al. 2012). Since statis- 
tical isotropy is predicted for the early Universe, analyses 
should focus on the primordial CMB, i.e. one from which 
secondary low-redshift cosmological signals have been re- 
moved. This should be done whether or not one believes 
the significance of the reported anomalies in the CMB. 

In this paper, we focus simultaneously on both 
the reconstruction method and the ISW effect as a 
means to estimate the primoridal CMB and test it for 
anomalies. We focus on three previously reported large- 
scale anomalies, namely: the low quadrupole power, the 
quadrupole/octopole alignment, and the octopole planarity. 

This work updates that of Francis & Peacock (2010) 
and Kim et al. (2012) in the following ways: 
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iLiKj Wo i^mpj 


Map 


Sparsely inpainting ILC W3 map 


TT C (\riTt\ 

llJKj vv O ^^liip ) 


iVicip 


Opdibciy iiipdiiiUiiig, lJ_jV_y wo iiidp 


Dela W5 finD) 


Map 


SiDPirsplv iriDJ^intinP' T^pl^^brnnillp AA/fS mpw) 


TT,C W7 Hnn^ 


Map 


^■p>Q tgpI'^;' in "r* q ittI" 1 n (T II ,( , \A/ / m q Tr\ 
Opdibcij' iiipdiiiUiiig, J.J_jV_y VV I iiidp 


TT C WQ (\-n'n\ 


A/rc»T^ 
iviap 


Opdibeiy iiipdiiiUiiig iijv_y VV c/ iiidp 


ZiViiT.OO ^^llip J 


iVldp 


Opdiociy iiipdiiiUiiig, ZiViiT.OO iiidp 


iy V OO l^liip ) 


iviap 


opdibeiy iiipdiiiUiiig iN V OO iiidp 


cmb _1 owl _a.lin _ i np a. i nt i ng 


code ylULi ) 


inpaints CMB or LSS maps (compatible Avith. ISAP^) 


zlVlAoo ioW products: 






io W -ZMAoo- i Uri 


Map 


TQAAT^ fT.-^vvi OA/r A QQ ^-^A TT^U 

ioW irom zMAoo and iUrl 


ID VV _ZiiViiT.DO_ VV ± 


iVidp 


T^AA/ from 9A/rA^^ anH AA71 

iO vv iiUiii ZiVliT.OO diiLl VV i 




iViap 


T^AA/ from 9A/rA^^ anri AA/'^ 

io VV irom zivii-voo ana wo 


TQAA7 9A/r A QQ AAT^i^ 
io VV _ZiVi/\oo_ VV 


Map 


TQAAT^ fi-om OA/I A QQ onrl AA^p; 

io VV irom ziViAoo and wo 


TQAA7 OA/r A QQ AA7i^T^/2>lo 

io W _zMAoo_vv oJJeia 


Map 


TQAA7 fr-orvi OA/I A QQ on/^ AA7c: T^-^lo 

ioW irom ZMAoo and Wo ueia 


TQAAT^ 9A/r A QQ AA77 
io VV -ZiViAoo- VV 1 


Map 


TQAA/^ fr-om OA/I A QQ on/^ AA7^7 

io VV irom ZiViAoo and vv ( 


TQAA7 OA/r A QQ AA7Q 
io W _zMAoo_VVy 


Map 


TQAAT^ fi-oTvi OA/I A QQ on/^ AA7Q 

ioW Irom ZiViAoo and wy 


ATA/'QQ TQAAT^ T->T.^^n/^+o. 

iNVoo ioW products: 






TQAAT^ MAAQQ 'TT^XJ 

io W _iN V oo_i Url 


Map 


TQAAT^ fT.-^vvi MAAQQ ^-^A TT^XJ 

io W Irom iN V OO and i Uil 


T^W NV^^ W1 

lO vv _iM V 00_VV ± 


A/Tar* 
iVi dp 


T^W from NV^^ p\nr\ W1 

lO vv iiUiii iM V OO diiLl VV ± 


io VV _iN V 00_ VV O 


iVidp 


T^AA/ from NA/^^ flnri AAl'^ 

io VV iiom iN V OO diia wo 


T^AA^ NV^^ AA^"^ 

iO vv _1N V 00_ VV O 


iVidp 


T^AA^ from NV^^ anH AA^^ 
iO VV iiOiii iN V OO diiQ wo 


ISW_NVSS_W5Dela 


Map 


ISW from NVSS and W5 Dela 


ISW_NVSS_W7 


Map 


ISW from NVSS and W7 


ISW_NVSS_W9 


Map 


ISW from NVSS and W9 


Statistics products: 






anomalies_1213 


code (F90) 


calculates quad/oct alignment and probability (requires HealPix) 


anomalies-octplan 


code (F90) 


calculates octopole planarity and probability (requires HealPix) 



Table 7. List of products made available in this paper in the spirit of reproducible research, available here: http: 
//www. cosmostat . org/anomaliesCMB .html. 



1. Tomography: We reconstruct the ISW signal due to both 
2MASS (Jarrett et al. 2000) and NVSS (Condon et al. 
1998) data (see Section 3) 

2. Trends: We look for trends by considering 11 different 
WMAP data sets (see Section 3), including 6 data sets 
for Avhich Ave reconstructed missing data using the sparse 
inpainting technique mentioned beloAV 

3. Sparse inpainting: Missing data in CMB and LSS maps 
are reconstructed using sparse inpainting (see Section 4, 
A08-SMF10), Avhich does not that assume the data are 
Gaussian or statistically isotropic. 

4. Tests for biases: We shoAv in Section 4 and Appendix 
A that the sparse inpainting method used is a bias-free 
reconstruction, Avhich does not introduce a spurious sta- 
tistical anisotropics nor a spurious ISW signal 

In this Avork Ave first investigate the impact of our bias- 
free sparse inpainting on various CMB maps and on the 
claimed CMB anomalies. We then subtract the ISW signal 
reconstructed in Section 4 and test again for impact on the 
three studied anomalies. Our main conclusions are made by 
identifying trends amongst the 11 CMB maps considered. 

Our conclusions are summarised in Table 6, namely: 

— The loAv quadrupole becomes more anomalous after 
sparse inpainting of CMB maps. 

— The quadrupole/octopole alignment becomes less 
anomalous after sparse inpainting of the CMB maps, 
contrary to Avhat Avas reported by Kim et al. (2012), 
Avho used a Gaussian constrained inpainting method. 



— Both the loAv quadrupole and the quadrupole/octopole 
alignment anomalies are reduced in significance after 
subtraction of the ISW signal, similarly to Avhat Avas 
reported by Francis & Peacock (2010). 

— We note that the reduced significance of the 
quadrupole/octopole alignment (both after sparse in- 
painting and ISW subtraction) is mainly due to changes 
in the quadrupole shape, not the octopole, suggesting 
that the main source of both anomalies could be the 
quadrupole. 

— We find that after inpainting and ISW subtraction, the 
octopole planarity is not significantly anomalous, con- 
trarily to the report of Francis & Peacock (2010) Avho 
found the octopole planarity had become more anoma- 
lous after ISW subtraction. We also note that the oc- 
topole planarity did not seem significant in any CMB 
map subsequent to the first year WMAP data. 

We therefore conclude that after application of sparse 
inpainting and subtraction of the ISW signal due to the 
loAv-redshift Universe out to 2; ~ 1, as estimated from the 
2MASS and NVSS surveys, that the Ioav quadrupole and the 
quadrupole/octopole alignment appear less significant and 
the octopole planarity remains unproblematic. In addition, 
it seems that both of these previously reported anomalies 
could be in fact due to the quadrupole only. Other hypothe- 
ses remain possible (e.g. exotic physics). 

In the spirit of participating in reproducible research, 
Ave make public all codes and resulting products Avhich 
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Map 


't' value 


Frobabiiity (/oj 


1 (Jrl 


0.9415 


9.60 


Wo 


U.9171 


1 /I o 

14.0 


Wo 


u.yiyu 


lo.o 


W7 


0.9261 


12.7 


ILC W9 


0.9o29 


11.8 




n SQ1 7 


/I 


Wo (mpj 


n ooon 

U.oozU 


ZD.l 


Wo (mpj 


0.8886 


24.4 


Wo Ueia (mpj 


o.8yiy 


Z0.4 


W7 (mpj 


0.9149 


15.1 


Wy (mpj 


0.9205 


lO.D 


rri/^TT TOTT7 

TOH -IbW 


0.8365 


38.7 


Wo -ibW 


0.8216 


42.8 


TT Tr' TOTT 7 

Wo -IbW 


r\ art A r\ 

0.8240 


42.0 


TT7^ TOTT7 

W7 -ibW 


0.8250 


41.8 


TT TTf/^V TOTT7 

ILC W9 -lb W 


0.8325 


39.8 


iKjri vv 1 (^mpj -low 


n 71 97 


7/1 1 
< 4.1 


vv O I liip ) -LkJ VV 


U. 1 uoo 


7(1 4 
1 u.^ 


W5 (inp) -ISW 


0.7115 


74.4 


W5 Dela (inp) -ISW 


0.7102 


74.8 


W7 (inp) -ISW 


0.7443 


65.4 


W9 (inp) -ISW 


0.7521 


63.4 


W3-ISW (due to 2MASS) 


0.9841 


1.6 



FPIO 

Table 5. The 't' value for the octopole as defined in de 
Oliveira-Costa et al. (2004) using nside=128^ calculated 
from the observed CMB maps {top) and after subtraction 
of the ISW field due to 2MASS and NVSS galaxies {mid- 
dle). The probability is determined from 1000 Monte-Carlo 
simulations and compared with results from FPIO (Francis 
& Peacock 2010, bottom row). 



Healpix software (Gorski et al. 2002; Gorski et al. 2005), ISAP^ soft- 
ware, the 2 Micron All-Sky Survey catalogue (2MASS) WMAP 
data the NRAO VLA Sky Survey (NVSS, Condon et al. 1998) 
and the Galaxy extinction maps of Schlegel et al. (1998). This work 
was supported by the European Research Council grant SparseAstro 
(ERC-228261) and by the Swiss National Science Foundation (SNSF). 



Appendix A: Testing statistical isotropy biases for 
sparse inpainting 

To justify that sparse inpainting is indeed an appropriate 
reconstruction method in the context of probing anomalies 
in the primordial CMB, we test whether we can indeed con- 
sider sparse inpainting as a bias-free reconstruction. Tto do 
this we test whether the three statistics studied in Section 
5 are unchanged after the application of sparse inpainting 
to CMB, LSS, and ISW simulated maps. 

Our aim here is to test whether statistically isotropic 
maps conserve this quality after sparse inpainting. It would 
also be interesting to test whether maps which are intrin- 
sically anomalous are still anomalous after sparse inpaint- 
ing. However, testing the latter requires a knowledge of the 
model creating the anomalies (exotic physics, foregrounds, 
etc. . . ) in order to test these realistically. In this appendix, 
we focus only on testing whether sparse inpainting creates 
spurious statistically anisotropic signatures in data where 
the true distribution is known to be statistically isotropic. 

We sparsely inpaint all CMB and LSS simulations using 
the public ISAP software^ of A08-SMF10 as described in 
section 4. 



Anomaly 


After Sparse 
Inpainting 


After ISW 
Subtraction 


Low Quad 


More Anomalous 


Less Anomalous 


Quad/Oct 


Less Anomalous 


Less Anomalous 


Alignment 






Oct planarity 


Less Anomalous 


Less Anomalous 



Table 6. Summary of results in this paper. We identify 
trends from the 11 WMAP data sets described in Table 
1 regarding how sparse inpainting and subtraction of the 
ISW signal due to 2MASS and NVSS galaxies can affect 
the three anomalies investigated. The third test for octopole 
planarity did not return significantly anomalous results in 
the first place. 



constitute the main results of this paper public. In Table 
7 we list the products which are made freely available 
as a result of this paper and which are available here: 
http : / / www . cosmostat . org/ anomaliesCMB . html. 
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referee for useful discussions. AR thanks Francois Lanusse, Florent 
Sureau, and Yves Revaz for computational help. We also thank 
Jacques Delabrouille for providing the wavelets ILC Syr map from the 
WMAP 5yr data. We used iCosmo ^ software (Refregier et al. 2011), 



http : //www . icosmo . org 



A.l. Impact of Sparse Inpainting on CMB data 

We start by considering the impact of sparse inpainting on 
CMB maps. We consider 1000 Gaussian random field reali- 
sations of a CMB map with an input power spectrum corre- 
sponding to the "vanilla" cosmology considered in Section 
5. We sparsely inpaint each simulation after applying the 
WMAP 7th year temperature analysis mask, as described 
in Section 3 and 4. 

We test the mean values for each statistic (low 
quadrupole, quadrupole/octopole alignment, and octopole 
planarity) before and after sparse inpainting. The results 
are reported in Table A.l and show that the sparse in- 
painting method introduces negligible biases in all three 
tests, compared the intrinsic cosmic variance. The test of 
quadrupole/octopole alignment can be somewhat altered 
due to sparse inpainting (see the standard deviation on the 
bias). 

Since with current data we observe a low quadrupole, 
we test whether sparse inpainting is also robust in cases 
where the quadrupole is low. We perform a new set of 400 
Gaussian random field simulations, where the input power 
spectrum is WMAP 7 data. Table A. 2 shows the results for 
the three statistical isotropy tests performed on this set of 
simulations before and after application of sparse inpaint- 
ing. Again, we see that biases are very small compared to 
cosmic variance. 



^ http : // j starck . free . f r/ isap . html 
http : //www . ipac . caltech . edu/2mass/ 
http://map.gsfc.nasa.gov 

http : / /heasarc . gsf c . nasa . gov/W3Browse/ all/ nvss . 
html 



14 



A. Rassat et al.: Removal of two large-scale CMB anomalies after subtraction of the ISW effect 



A.2. Impact of Sparse Inpainting on LSS data 

For the reconstruction of the ISW signal, we must first re- 
construct full-sky LSS maps (Section 4). Thus we also want 
to test whether the statistics of the LSS maps are altered 
by sparse inpainting. To do this, we specifically test how 
Gaussian random field simulations of 2MASS and NVSS- 
like surveys are affected by sparse inpainting using the same 
masks as used in Section 3. The simulated maps are for the 
galaxy overdensity (i.e. including the galaxy bias). Results 
are shown in Tables A.3 (2MASS) and A.4 (NVSS). Again, 
biases are negligible compared to cosmic variance for all 
three statistics considered. 

A.3. Impact of Sparse Inpainting on the ISW signal 

Finally, we want to test whether sparse inpainting can alter 
the reconstructed ISW signal, either by modifying tests of 
statistical isotropy or by introducing a spurious ISW signal. 
We consider 1000 Gaussian random realisations of NVSS 
galaxies and 1000 uncorrelated WMAP 7 realisations. The 
goal is to see if inpainting of each map with its respec- 
tive masks can introduce a spurious ISW signal, when none 
should be measured. The results are shown in the first line 
of Table A. 5, which gives the value of the cross-correlation 
quadrupole (C^t/=2)- Since the maps are uncorrelated by 
construction, we expect a median cross-correlation equal to 
zero, which is what we find both before and after sparse in- 
painting. This indicates that our method does not introduce 
a spurious cross-correlation signal. 

We also test the quadrupole/octopole alignment and 
octopole planarity of the reconstructed ISW signal from 
NVSS galaxies (where the ISW map is constructed from 
NVSS and CMB simulations with different masks). We find 
that the tests are unbiased after sparse inpainting. 





ILC Th 
simulations 

a 


ILC Th 
inpainted 
simulations 
b 


a — b 


Expected 
Quadrupole 


1269.0 ±806.0 


1229.9 ± 855.1 


39.1 ±334.3 


Expected 
(n2 • fis) 


0.489 ± 0.302 


0.503 ± 0.296 


-0.014 ±0.181 


Expected 
(t) value 


0.78 ±0.117 


0.78 ±0.119 


-0.0012 ±0.061 



Table A.l. Testing whether sparse inpainting (A08- 
SMFIO) affects the tests of statistical isotropy used in 
Section 5 for 1000 simulated Gaussian random fields with 
input power spectrum taken as the best-fit theoretical 
power spectrum for WMAP7 (i.e. with no low quadrupole). 
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simulations 
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ILC W7 
inpainted 
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a — b 


Expected 
Quadrupole 


OKO 1 _L 1 n 
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Table A. 2. Testing whether sparse inpainting (A08- 
SMFIO) affects the three tests of statistical isotropy used 
in Section 5 for 400 simulated Gaussian random fields with 
input power spectrum taken as WMAP7-like data (i.e. with 
low quadrupole). 





2MASS gal 
simulations 

a 


2MASS gal 
inpainted 
simulations 
b 


a — b 


Expected 

Quadrupole 

[xlO-^] 


2.96 ± 1.93 


3.00 ±2.49 


-0.041 ± 1.44 


Expected 
(n2 • hs) 


0.546 ± 0.290 


0.536 ±0.303 


0.011 ±0.392 


Expected 
(t) value 


0.782 ±0.120 


0.801 ±0.115 


-0.019 ±0.145 



Table A.3. Testing whether sparse inpainting (A08- 
SMFIO) affects the tests of statistical isotropy used in 
Section 5 for 1000 simulated Gaussian random fields with 
input power spectrum taken as the theoretical galaxy power 
spectrum for 2MASS galaxies. 
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